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Summary 

Broadband shortwave and longwave radiance 
measurements obtained from the Nimbus 7 Earth 
Radiation Budget scanner were used to develop re- 
flectance and emittance models for the Sahara, Gib- 
son, and Saudi Deserts. These deserts will serve as 
in-flight validation targets for the Earth Radiation 
Budget Experiment being flown on the Earth Radia- 
tion Budget Satellite and two National Oceanic and 
Atmospheric Administration polar satellites. The di- 
rectional reflectance model derived for the deserts 
was a function of the sum and product of the cosines 
of the solar and viewing zenith angles, and thus reci- 
procity existed between these zenith angles. The 
emittance model was related by a power law of the 
cosine of the viewing zenith angle. 

Introduction 

Three Earth Radiation Budget Experiment 
(ERBE) missions will be launched during 1984-1986 
to study the regional, zonal, and global radiation 
properties of the Earth. To meet the stringent ac- 
curacy goals set for ERBE, the satellite orbits must 
permit frequent coverage of the entire globe, and the 
instruments must provide high radiometric accuracy. 
Coverage errors are basically fixed by the number of 
satellites and their orbits (Harrison et al. 1983), while 
radiometric errors are determined by the design of 
the sensors and the degree to which their calibra- 
tions are known in flight. All the ERBE instruments 
were extensively calibrated and characterized in a 
ground test facility; nonetheless, onboard calibration 
systems were provided to recalibrate the sensors and 
to validate their long-term stability in flight. A more 
complete description of the ERBE sensors and their 
ground and onboard calibration systems is given by 
Barkstrom (1984). 

The need for in-flight calibration became appar- 
ent in the early 1960’s when radiometers exposed 
to the harsh space environment often degraded at 
alarming rates and in unknown manners. Conse- 
quently, onboard calibration systems were added that 
generally took the form of solar diffusers, incan- 
descent lamps, or both for shortwave sensors and 
temperature-controlled blackbodies for longwave sen- 
sors (Williamson 1977). Unfortunately, in-flight cali- 
bration systems can also degrade, fail, or become un- 
usable for numerous reasons. Futhermore, their ap- 
plicability is questionable because they are usually 
activated in special calibration sequences that are 
substantially different from their operational Earth- 
viewing modes. 

One in-flight calibration technique that eliminates 
some of the problems encountered with onboard cal- 


ibration systems is the use of operational measure- 
ments that view natural Earth targets as calibration 
values. Selected Earth sites, in effect, become the 
solar diffuser plates and thermal blackbodies. For 
instance, the bright gypsum sands and the nearby 
dark lava flows at White Sands, New Mexico, are 
occasionally used as relative comparison targets for 
satellite sensors. The primary use of such sites for 
the ERBE missions will be to serve as a long-term 
stable target for calibrations and for comparisons 
of measurements (1) among the three ERBE mis- 
sions, (2) between different ERBE instruments on a 
given mission, (3) between ERBE results and those 
of previous missions (Nimbus 2 and 3, medium res- 
olution infrared radiometer; Nimbus 6 and 7, Earth 
Radiation Budget (ERB) scanner), and (4) between 
ERBE results and those of future radiation budget 
missions. 

The present paper is concerned with, first, the se- 
lection of target sites amenable to the ERBE mis- 
sions, and then, the directional modeling of the 
broadband reflection and emission properties of those 
sites. For a number of reasons given later, all the se- 
lected sites were deserts: Sahara, Gibson, and Saudi. 
Nimbus 7 ERB scanner data were used to develop 
the cloud-free desert models, and comparisons are 
made between the Nimbus 6 and 7 ERB results for 
the Saudi Desert. 


Symbols 


A albedo 

ClWiCsw longwave and shortwave azimuthal 
phase function constants 

E^w longwave radiant exitance, W-m -2 

Flw^Fsw measured longwave and shortwave 
radiant fluxes, W-m -2 

LlWi Fsw broadband longwave and shortwave 
radiance measurements, W-m~ 2 -sr -1 

Liw directional emittance, W-m~ 2 -sr“ 1 

(see eq. (9)) 

LlW6 directional emittance to zenith angle 
0, W-m _2 -sr _1 

M longwave regression coefficient (see 

eq. (10)) 

N shortwave regression coefficient (see 

eq. (3)) 

Plw -> Fsw longwave and shortwave azimuthal 
phase functions 



R 

bidirectional reflectance (see eq. (1)) 

R 

directional reflectance (see eq. (2)) 

^-min 

minimum bidirectional reflectance 

r a ,r s 

Earth radii to top of atmosphere and 
to satellite 

S 

Earth-Sun distance-corrected solar 
flux, W-m -2 

S o 

mean Earth-Sun distance solar flux, 
1376 W-m" 2 

U,Uo 

cosines of viewing and solar zenith 
angles, respectively 

U max 

cosine of maximum viewing zenith 
angle observed by MFOV sensor 

V,V 0 

sines of viewing and solar zenith 
angles, respectively 

X 

UUq/{U + Uo) or (sec0 + sec 0 q) _1 

Y 

RUUo or RUUq/Psw 

YqM 

shortwave regression coefficients (see 
eq. (3)) 

a 

cone angle of MFOV sensor 

a ma.x 

maximum cone angle of MFOV sensor 

0 

clock angle of MFOV sensor 

7 

scattering angle (see eq. (4)) 

9, 6q 

viewing zenith angle and solar zenith 
angle, respectively 

a 

standard error of estimate for Y 


solar azimuth angle 

Abbreviations: 

AVHRR 

advanced very high resolution 
radiometer 

ERB 

Earth Radiation Budget (sensors 
onboard Nimbus 6 and 7) 

ERBE 

Earth Radiation Budget Experiment 

ERBS 

Earth Radiation Budget Satellite 

HIRS 

high resolution infrared radiation 
sounder 

IFOV 

instantaneous field of view 

LST 

local standard time 

MFOV 

medium field of view 


NFOV narrow field of view 

NO A A National Oceanic and Atmospheric 

F,G Administration Tiros N-series 

satellites 

WFOV wide field of view 

Site Selection 

Required and Desired Site Properties 

The rationale for the selection of the present 
target sites is given below. Certain site properties 
must be classified as requirements, while others are 
merely desirable properties. 

Stability. The most basic requirement of any tar- 
get site is that it must provide a stable, predictable 
radiometric signal when observed from space. A site 
must offer long-term stability with no appreciable 
surface changes expected during the period of in- 
terest in order for direct comparisons to be made 
between previous as well as future radiation budget 
measurements and those for ERBE. Even short-term 
variances due to ice, snow, wet soil, and vegetation 
growth, for instance, are undesirable features in that 
ERBE target sites could be utilized on a daily basis, 
when cloud-free. Although all areas of the globe ex- 
perience some degree of cloudiness, strong weighting 
was given to those sites where cloudiness and precip- 
itation are minimal. 

High albedo/low solar zenith angle. All the ERBE 
shortwave sensors view space and/or the dark side of 
the Earth during flight and thus experience a zero 
input condition. The most appropriate additional 
reference condition would be one near the maximum 
levels measured during flight, which occur for high- 
albedo sites at low solar zenith angles. Clouds, 
ice, and snow have high albedos but do not provide 
stable, predictable radiometric surfaces. In general, 
deserts are stable and predictable and have high 
albedos, and several large deserts are located at low 
latitudes where the lowest possible solar zenith angles 
will be experienced by each of the ERBE satellites. 

Size . It would be desirable to have target sites for 
all three of the ERBE sensors, which are the narrow- 
field-of-view (NFOV) scanner and the medium-field- 
of-view (MFOV) and wide-field-of-view (WFOV) 
nonscanners. However, no acceptable site is suffi- 
ciently large for the WFOV sensor, which views the 
entire Earth disc, and its in-flight calibration will 
rely heavily on the onboard calibration systems. The 
MFOV sensor views an Earth-central-angle nadir 
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scene of 10°, and for observations each day, the site 
must be about 10° in latitude by 36° in longitude. 
The NFOV scanner has an IFOV of about 3° by 3°, 
and its target sites could be relatively small. How- 
ever, scanner measurements will be assembled by re- 
gions that are 2.5° in latitude and longitude, which 
thus become convenient site dimensions. Approxi- 
mately 20 (limb) to 120 (nadir) NFOV scanner mea- 
surements will be obtained from regional sites during 
each day and night overpass. 

Uniformity . There is no requirement that a tar- 
get site be radiometrically uniform provided that the 
entire site is sampled during each orbital passage. 
However, uniformity is certainly a desirable site prop- 
erty in that histograms of the shortwave and long- 
wave radiances (from AVHRR, HIRS, and/or ERBE 
NFOV) will be used to detect the presence of clouds 
and, perhaps, wet sand and dust storms. Lambertian 
or near-Lambertian target sites are desirable in that 
they can be more accurately modeled. 

Potential Areas 

Only low-latitude deserts meet all the require- 
ments: long-term stability and predictability, high 
albedos combined with low solar zenith angles, and 
large size. Furthermore, most of these deserts are lo- 
cated in northern Africa, the Arabian Peninsula, and 
central Australia, all of which have low-cloudiness 
and low-precipitation statistics (Sherr et al. 1968). 
Nimbus 6 ERB scanner data taken during June 1977 
(scanner pointed at nadir) were used to search these 
large areas for potential target sites. In particular, 
the search was directed toward finding near-uniform 
regions that might best serve as NFOV target sites. 

Selected Target Sites 

The area selected as the MFOV target site is lo- 
cated in the Sahara Desert and is bounded by lat- 
itudes 17°N to 27°N and longitudes 3°W to 33°E. 
This site can generally be characterized as having 
vast expanses of sand dunes, sand seas, arid gravel, 
and rock outcroppings with very little moisture or 
vegetation (Short et al. 1976). Nadir radiances from 
Nimbus 6 revealed maximum variances of ±20 per- 
cent for the entire site, but maximum variances of 
only ±10 percent for more than 90 percent of the 
site. Because of solar zenith angle problems, ex- 
plained later, a second MFOV site was considered 
which spanned the entire central Australian Con- 
tinent. This area has large variances (±30 per- 
cent) plus seasonal vegetation growth (mostly range 
grasses) and was not considered acceptable as a tar- 
get site. 


The region selected as the primary NFOV target 
site is located in the Ar Rab al Khali Basin of 
the Arabian Peninsula and is bounded by latitudes 
20.0°N to 22.5°N and longitudes 50.0°E to 52.5°E. 
This general area is often referred to as the Arabian 
Empty Quarter, but because it is located totally 
within the borders of Saudi Arabia, it will be referred 
to here as the Saudi Desert. This site is almost 
entirely sand dunes and sand seas, and the lack of 
moisture and the saline nature of the sand cause it to 
be virtually free of all vegetation (Short et al. 1976). 

For two opportunities during each 72-day preces- 
sion cycle of the Earth Radiation Budget Satellite 
(ERBS), the minimum solar zenith angles for the 
satellite as it passes over the Saudi and Sahara sites 
will range from 0° to 45° depending on the time of 
year (i.e., solar declination). Solar zenith angles will 
range from 35° to 58° for the NOAA F (1430 LST) 
and from 55° to 74° for the NOAA G (0800 LST) 
Sun-synchronous spacecraft as they pass over these 
sites during the course of a year. While the ranges 
of solar zenith angles for the ERBS and NOAA F 
spacecraft are acceptable, the maximum angles for 
the NOAA G spacecraft, which occur during the fall 
and winter seasons at the Sahara and Saudi sites, are 
somewhat large and will not provide the high radio- 
metric values desired for shortwave calibrations. This 
problem can be alleviated by the addition of sites 
located in the Southern Hemisphere for use during 
October through March. 

The region selected as an alternate NFOV target 
site is located in the Gibson Desert of western Aus- 
tralia and is bounded by latitudes 25. CPS to 27.5°S 
and longitudes 120. 0°E to 122. 5°E. This site has a 
mixture of sand dunes, dry salt lakes, rock outcrop- 
pings, and arid steppe vegetation (Short et al. 1976). 
However, after infrequent heavy rains, the lakes tem- 
porarily flood, and the sparse vegetation flourishes. 
These occurrences are undesirable from the present 
standpoint because they would prevent use of the site 
until arid conditions returned. 

Broadband shortwave, Lg\y and longwave, 
radiance measurements from the Nimbus 6 ERB 
scanner are presented in figure 1 to illustrate the 
surface uniformities of the Saudi and Gibson NFOV 
target sites. These nadir measurements represent 
100-km swaths that passed directly across the centers 
of the sites as the spacecraft coursed northward 
on daylight passages. The shortwave and longwave 
radiance variances for these as well as the adjacent 
swaths (i.e., total coverage) were less than 5 percent 
for both NFOV sites. 
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DEGREES LATITUDE NORTH OF SOUTHERN BOUNDARY 

Figure 1. Illustration of uniformity for NFOV scanner sites. 

Data from Nimbus 6 ERB scanner; 0 — 6°; June 1977. 

Satellite Measurements 
Satellites and Instruments 

Radiance measurements from the Nimbus 7 ERB 
scanner were used for modeling purposes, and they 
were compared with the Nimbus 6 results for the 
Saudi Desert. The Nimbus 6 satellite was launched 
June 12, 1975, to an altitude of 1100 km and an 
inclination of 100°, and the Nimbus 7 satellite was 
launched October 24, 1978, to an altitude of 950 km 
and an inclination of 99° — both into noon, Sun- 
synchronous orbits. 

The ERB scanner had four telescopes, each of 
which had a broadband shortwave (0.2-4 /im) chan- 
nel and a broadband longwave (5-50 /im) channel. 
It scanned from horizon to horizon along the or- 
bital groundtrack and to a zenith angle of 72° in the 
crosstrack direction. This multiaxis scanning capa- 
bility provided measurements over a wide range of 
angular parameters which are particularly useful for 
directional modeling purposes. At the nadir, the spa- 
tial resolution (Nimbus 7) was about 90 km, and it 
increased to about 250 km at the horizon. A more 
detailed description of the ERB instrument and its 
operation is given by Smith et al. (1977). 

Data Sets 


longwave radiances from the Sahara Desert measured 
by the Nimbus 7 ERB scanner from November 1978 
to May 1980 (Stowe and Fromm 1983). These data 
were sorted into 191 ranges of angles, called bins, to 
provide discrete parameter values for modeling and 
display. Bins were sorted by viewing zenith angle, 
0, azimuth angle, (f> (<j) = 0° for forward reflection; 
<j) = 180° for backward reflection), and E/o, cosine of 
the solar zenith angle. Bin angular ranges are given 
in table I. For each of the data bins, the bin de- 
scription (#, </>, E/q), the mean of both Lg\y/Uo and 
Llw, and the bin population were available in tab- 
ular form. In addition, the albedo, A , was computed 
for each of the four E/o ranges by using a mean Earth- 
Sun distance solar flux, So, of 1376 W/m 2 . This data 
set was compiled and processed under the direction of 
V. R. Taylor and L. L. Stowe of the National Oceanic 
and Atmospheric Administration/National Environ- 
mental Satellite, Data, and Information Service. 

The Nimbus 7 data used to model the Gibson 
and Saudi Deserts were taken during several days in 
December 1978 and March, July, and December 1979 
when the spacecraft passed near the sites. These 
measurements were individual shortwave and long- 
wave radiance pairs obtained from within the bound- 
aries of the 2.5° regional sites. The December 1978 
through July 1979 data provided maximum, medium, 
and minimum solar zenith angles at both sites, and 
the December 1979 data were added primarily for 
direct comparisons with the December 1978 data to 
determine whether sensor degradations had occurred 
during the year. Nimbus 6 data for the Saudi Desert 
taken during August 1975 (Suttles 1981) were in- 
cluded for comparison with the Nimbus 7 results. 
Data set populations and angular ranges are given 
in table II. 

Measurements taken with clouds within the in- 
stantaneous fields of view of the ERB scanner were 
eliminated from the Sahara data set by using the win- 
dow radiances (11 /im) of the temperature-humidity 
infrared radiometer (THIR, also onboard Nimbus 7) 
in an editing technique developed by Stowe (1983). 
Cloudy days were eliminated from the Gibson and 
Saudi data sets by analyzing the general levels of the 
ERB Lsw an d Li\y values (overcast) and by re- 
quiring continuity in the daily plots of L^\y versus 0 
(partial cloudiness). 

Reflectance Model 

Sahara Directional Model 

Bidirectional reflectance is defined here as 


The data set used to develop the basic desert 
models consisted of 123 141 pairs of shortwave and 


R = 


T^SW_ 

SUo 


( 1 ) 
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Figure 2. Correlation of directional reflectance for Sahara 
Desert as a function of viewing and solar zenith angles. 

and is a function of 9q, and <f>. ( S is the distance- 
corrected solar flux.) Azimuthal dependence is elimi- 
nated for the directional reflectance, which is defined 
as 

— If 77 

R=- Rd(f> (2) 

n Jo 

Barkstrom (1973), using single-scattering approxi- 
mations, found that RUUq could be expressed as a 
function of UUq/(U + Uq) where U = cos 9. Staylor 
(1985) also found these parameters to be useful for 
the directional modeling of overcast low-, middle-, 
and high-water clouds and high-ice clouds. 

Directional reflectance parameters, RUUq, for the 
Sahara Desert data are presented in figure 2 as a 
function of UUq/(U + Uq), which is the inverse of 
the total atmospheric slant path, sec 9 + sec0Q. The 
data are well correlated with a curve related by the 
equation 

Y = Y 0 + Y,X n (3) 

where X = UU 0 /(U + U Q ) and Y = RUUq. A 
general discussion of the technical significance of the 
regression coefficients Yo> Fl, and N is given by 
Staylor (1985). Briefly, Yq accounts for the nonzero 
radiance due to atmospheric refraction and scattering 
measured near the terminator (Uo = 0, X = 0). The 
coefficient Y\ relates the general brightness level of 
the surface: 0.1 to 0.2 for clear oceans, 0.2 to 0.5 for 
soils, 0.6 to 1.1 for deserts, 1.2 to 1.5 for clouds, and 
1.8 to 2.4 for snow. The coefficient N has several 
components, one of which relates the transmission 
within the surface medium, primarily grains of sand 


in the present case, and another of which relates the 
transmission of the atmosphere above the surface. 

Azimuthal Phase Function 

Directional reflectance is a useful parameter for 
correlative purposes and for the calculation of albedo. 
However, during an overflight of a target site, the 
ERBE NFOV scanner will not observe the site at 
sufficient azimuth angles for the accurate computa- 
tion of R. Therefore, directional reflectance has little 
direct application to the present task of estimating 
target radiances unless the site reflectance is isotropic 
(i.e., R = R). It is clear that an azimuthal phase 
function is needed that will convert the modeled R 
values into the required R values. The azimuthally 
binned Sahara data set was ideally suited to this task. 

Staylor (1985) found that the bidirectional re- 
flectances for clouds were well correlated with scat- 
tering angle, 7, which can be expressed as 

cos 7 = sin 9 sin 9q cos (f) — cos 9 cos 9q (4) 

where 7 = 0° is forward scattering, and 7 = 180° is 
backward scattering. Plots of R versus 7 for constant 
bin values of 6 and 9q indicated that the Sahara 
reflectances behaved in a Rayleigh-like manner and 
could be approximated as 

R ~ -Rmin(l + Cs\V C °S 2 l) (5) 

where iZmin is a function of 0 and 9q, and C$w ls a 
constant. The azimuthal phase function, Psw, can 
be expressed as 

= R = 1 + CswjUUp — Wq cos 0 ) 2 
R l + C SW [(UUo) 2 + l/2(VV 0 ^} [> 

where V = sin 0 and Vq = sin 9q. 

Equation (6) indicates that azimuthal anisotropy 
should peak when the VVq term is at its maximum 
value. Measured values of R/R are plotted in figure 3 
as a function of 7 for six U, Uq bins which have 
large VVq values (i.e., small values of UUq). These 
data are well correlated with a curve generated by 
equation (6) with mean values of U and Uq, and 
Csw ~ 0.33. 

The azimuthal phase function, as related by equa- 
tion (6), has several notable features. First, it will ac- 
commodate predominately backward-reflecting sur- 
faces ( Csw > 0) such as the Sahara Desert, predom- 
inately forward- reflecting surfaces (Cgw < 0), and 
isotropic surfaces (Csw = 0)- Second, azimuthal de- 
pendence is minimal if one or both of the zenith an- 
gles are near their nadirs (VVq « 0), which is a gen- 
erally observed phenomenon. And third, the phase 
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Figure 3. Illustration of Rayleigh-like shape of azimuthal 
phase function as a function of scattering angle. Data 
from Sahara Desert for high zenith angles. 


function exhibits zenith angle reciprocity, which is 
discussed in more detail later. 

Gibson and Saudi Directional Models 

Bidirectional parameters RUUq for the Gibson 
Desert are presented in figure 4 (top, left scale) as 
a function of X. These data were taken during 
July 1979, when the noontime solar zenith angle 
was near its maximum yearly value (6q « 50°), and 
best illustrate the relatively high anisotropy of this 
target site. The backward-reflecting data lie along 
the top dashed line, the forward-reflecting data lie 
along the bottom dashed line, and the side-reflecting 
data lie between. However, when normalized by 
the azimuthal phase function, PswiPsw — 0.60), 
these data collapse into a single line (bottom, right 
scale). All the Gibson data for the 4 sampled 
months are given in figure 5 and are well correlated 
with a regression curve related by equation (3). 

Nimbus 7 .Saudi Desert data for the 4 sampled 
months are presented in figure 6 (top, left scale) as 
a function of X. Similarly, the Nimbus 6 Saudi data 
for August 1975 are also given (bottom, right scale). 
Both data sets were normalized by Pswi^SW — 
0.18) and show good agreement with their regression 
curves. All regression coefficients used for the desert 
reflection models are listed in table III. 



RUIL 


sw 


Figure 4. Normalization of bidirectional reflectance for 
Gibson Desert site, 0q m 50°, July 1979. Note difference 
in left and right vertical scales. 


Anisotropic reflection diagrams (i.e., R/A) are 
given in figure 7 for the Gibson (left) and Saudi 
(right) Deserts. The diagrams show that both desert 
sites are predominately back-reflecting, and that the 
Saudi Desert has the more desired, near-isotropic 
pattern. 

Albedo 

Local albedo can be expressed as 

A — 2 f 1 RUdU (7) 

Jo 

which with the substitution of equation (3) becomes 

-4 = 2 W + du (8) 

The albedo curves given in figure 8 were computed 
from equation (8) with the Nimbus 7 coefficients 
listed in table III. Bin-integrated albedos for the Sa- 
hara Desert (Taylor et al. 1983) are in good agree- 
ment with their computed curve. The Saudi site is 
clearly the brightest of the three deserts, followed by 
the Sahara, with the Gibson site being the darkest 
desert. 
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Figure 5. Correlation of directional reflectance for Gibson Figure 6. Correlation of Nimbus 7 (left scale) and Nimbus 6 
Desert site as a function of viewing and solar zenith (right scale) directional reflectances for Saudi Desert 

angles. site. 



Figure 7. Anisotropic reflectance diagrams (R/A) computed Figure 8. Albedos computed from Saudi, Sahara, and Gib- 

from Gibson and Saudi Desert models. Uo = 0.6. son Desert models. Symbols are bin-integrated values 

for Sahara Desert from Taylor et al. (1983). 
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Figure 9. Correlation of limb-darkening with zenith angle 
for Sahara Desert. 

Emittance Model 

There was little variation of the emitted longwave 
radiances with azimuth angle for the Sahara data set, 
and the directional emittance defined as 

1 f ^ 

L lw — ~ Llxy d(j) (9) 

7T Jo 


was computed for all Mo bins. Values of L lw 
plotted in figure 9 versus sec 6 on a log-log scale are 
quite linear for each value of Uo and can be related 
as a power law of U as 

L lw,0 ~ P lw, o°i sec @) ^ = Llw^oU ( 10 ) 

Staylor (1985) previously found that this limb- 
darkening function also applied for overcast clouds. 

Initial attempts to model the Plw data for the 
Gibson and Saudi Deserts in a similar manner were 
somewhat disappointing. For instance, the Gibson 
longwave radiances for July 1979 ($o ~ 50°) behaved, 
on a reduced scale, in a manner very similar to the 
shortwave reflectances shown in figure 4. Namely, 
the backward radiances (</> = 135° to 225°) were 
higher and the forward radiances ( <j> = 315° to 45°) 
were lower than the mean values. It was found 
that this emission anisotropy could be removed, to a 
large degree, by an azimuthal phase function, shown 
below, that was identical in form to that used for the 
shortwave anisotropy 


Plw = 


l LW 

Plw 


1 + C LW {UU 0 - Wo cos 4>) 2 
l + C LW [(UUo) 2 + l/2(VVo) 2 } 

( 11 ) 



Figure 10. Correlation of limb-darkening with zenith angle 
for Gibson and Saudi Desert sites. 

The normalized longwave radiances, Plw/ p LW j 
for the Gibson and Saudi sites, presented in figure 10, 
are generally linear even to quite large zenith angles. 
All regression coefficients and constants determined 
for the desert emission models are listed in table IV. 
The longwave radiant exitance expressed as 

O'jr _ 

e lw = 2 + ^ l w,o° ( 12 ) 

can be computed from these values. 

Discussion of Results 

Reflectance Model 

Comparisons of the December 1978 and Decem- 
ber 1979 Nimbus 7 data for both the Gibson and 
Saudi Deserts, given in figures 5 and 6, respectively, 
indicate that there were no significant differences in 
the levels of the measurements. Therefore, it was 
assumed that the sensor responses were constant 
during this period and during the somewhat longer 
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period of the Sahara data set. Stability is impor- 
tant, particularly for modeling purposes, in that a 
continuous degradation of the sensors during the De- 
cember 1978 to July 1979 period, for instance, could 
have been erroneously incorporated into the models 
as an effect of 6q. 

The multiaxis ERB scanner obtained measure- 
ments for all azimuth angles; however, with the as- 
sumption of symmetry about the principal plane, the 
data for 0 = 0° to 360° were mapped into 0 = 0° 
to 180° for the binned Sahara data set. Desert sand 
dunes are usually aligned with the prevailing winds, 
and this assumption might not be valid for such sur- 
faces. Although the validity of the symmetry as- 
sumption could not be assessed for the Sahara data 
set, there were no significant differences between the 
measurements for (j) = 0° to 180° and (j) = 180° to 
360° for the Gibson and Saudi Deserts. 

The regression parameters Y$/Y\,N, and Csw 
for the Saudi Desert (table III) were essentially iden- 
tical for the Nimbus 6 and 7 ERB scanners, but the 
Y\ value for the Nimbus 6 was about 9 percent higher 
than that for Nimbus 7. The inability to determine 
which of these two values is correct or more nearly 
correct emphasizes the fact that target sites can only 
be used to assess long-term sensor stability and not 
absolute sensor accuracy. 

Error sources for the present modeled data in- 
clude site nonuniformity, instrument noise, me- 
teorological variations, and modeling errors. In 
spite of such potential errors, the model dispersions 
(Standard error of the estimate/ Mean) were only 
5 percent to 8 percent (table III), and site nonunifor- 
mities such as those shown in figure 1 could account 
for a major portion of these dispersions. 

All zenith angle terms in the directional model, 
equation (3), and in the azimuthal phase function, 
equation (6), appear only in sum or product com- 
binations (i.e., U + Uo.UUo, and VVo), and there- 
fore, reciprocity exists between the viewing and so- 
lar zenith angles. (The values of 0 and do can be 
switched, and the value of R will be unchanged.) 
This feature of the target sites allows the present lim- 
ited range of Oo values to be extended to the more 
extensive present range of 6 values. This is a use- 
ful feature in that 6q will often exceed 50° during 
passes of the satellite over the sites, and it also lends 
a degree of confidence to the extrapolated values of 
the site albedos given in figure 8. The general form 
of the present model is probably applicable to most 
deserts and arid regions, and perhaps to other diffuse 
surfaces. 

The stated reason for including the Gibson Desert 
as an NFOV site was to provide a high-radiance tar- 
get for the NO A A G scanner from October to March 


when the solar flux impinging on the surface, SUo , 
at the Saudi site would be low. Unfortunately, the 
albedo of the Gibson Desert is considerably less than 
that of the Saudi Desert (fig. 8), meaning that its 
shortwave radiances would exceed those for the Saudi 
Desert only during December and January. This cal- 
culation is based on a nominal 0800 LST equatorial- 
crossing time, but should the NO A A G spacecraft be 
launched and/or precess to earlier crossing times, use 
of the Gibson site would become a necessity. 

Unlike the NFOV sensors, the MFOV sensors do 
not measure radiance, per se, but instead measure 
radiant flux, F. The shortwave flux can be related in 
terms of the shortwave radiance as 

C 2tt rot max 

Fsw — / / Lsw s ^ n a cos &dad(3 (13) 

Jo Jo 

where a and /3 are the polar angles of the MFOV 
sensor. With substitution, this expression becomes 


Fsw 





dUd(3 

(14) 


where r a and r s are the earth radii to the top of 
the atmosphere and to the satellite, respectively. 
Equation (14) can be numerically integrated by using 
the scanner-derived regression coefficients to provide 
an estimate of Fg\y for the Sahara site. 


Emittance Model 

The longwave radiance model given by equa- 
tions (10) and (11) provided good estimates for the 
noontime Nimbus 6 and 7 ERB measurements, but 
its regression values, L^wfi 0 and M apply only to 
noontime and not to other hours of the day. For in- 
stance, the noontime values of Llw,0° and E^w are 
quite linear with Uq for the Sahara Desert (table IV), 
yet these values are larger in the afternoon than in 
the morning for the same Uq. The precessing ERBS 
scanner should provide sufficient longwave measure- 
ments to model the diurnal and seasonal variations 
of L^iy o 0 and M. 

The azimuthal anisotropy for the longwave ra- 
diances came somewhat as a surprise. However, 
it should be expected that the side of sand dunes, 
rocks, and other vertical protuberances facing to- 
ward the Sun would be somewhat warmer and thus 
would emit more energy in that direction than would 
the opposite side in its direction. These protuber- 
ances may also be the major cause of the short- 
wave anisotropy. (Note that the Cjjyy values given 
in table IV are an order of magnitude smaller 
than the Csw values given in table III.) Longwave 
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anisotropy can be avoided by taking radiance mea- 
surements either during nighttime orbital passes af- 
ter the thermal anisotropy has dissipated or at the 
nadir (V = 0 ,Plw = 1) during daytime orbital 
passes. 

The MFOV longwave flux for Plw = 1 can be 
expressed as 

= (15) 

At intervals during the course of the ERBE mission, 
NFOV scanner measurements of the Sahara site that 
lie within the fields of view of the MFOV sensors 
will be modeled to provide the regression coefficients 
needed for equations (14) and (15). This procedure 
will provide both shortwave and longwave compar- 
isons between measurements of the NFOV scanner 
and MFOV sensors. Such comparisons were not pos- 
sible for ERB because the instrument packages did 
not have MFOV sensors. 


Concluding Remarks 

Broadband shortwave and longwave radiance 
measurements obtained from the Nimbus 7 Earth 
Radiation Budget (ERB) scanner were used to de- 
velop cloud-free reflectance and emittance models 
for the Sahara, Gibson, and Saudi Deserts. These 
deserts will serve as in-flight validation targets which, 
together with onboard calibration systems, will be 
used to provide radiometric stability for the Earth 
Radiation Budget Experiment (ERBE) sensors dur- 
ing the course of their mission lifetimes. The large 
Sahara Desert will serve as the target site for the 
medium-field-of-view (MFOV) nonscanner sensors. 
The smaller Saudi Desert will serve as the primary 
target site for the narrow-field-of-view (NFOV) scan- 
ner sensors, with the Gibson Desert as the alternate 
NFOV site. These deserts were selected as target 
sites because of their stability, clear skies, high albe- 
dos, low latitudes (i.e., low solar zenith angles), size, 
and uniformity. 

The directional reflectance model for these desert 
sites was a function of sum and product terms of 
the cosines of the solar and viewing zenith angles, 
and thus reciprocity existed between these zenith an- 
gles. This feature of the sites allows the present lim- 
ited range of solar zenith angles to be modeled to 
the more extensive present range of viewing zenith 
angles. The emittance model was related by a 
power law of the cosine of the viewing zenith an- 
gle. The general forms of both the reflectance and 
emittance models are probably applicable to most 


deserts and arid regions, and perhaps to other dif- 
fuse surfaces. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
December 5, 1985 
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TABLE I. SAHARA DATA BINS 


Variable and unit 

Data bin ranges 

9, deg 

“6-15 


15-27 


27-39 


39-51 


51-63 


63-75 


75-87 

<f>, b deg 

0-9 


9-30 


30-60 


60-90 


90-120 


120-150 


150-171 


171-180 

Uq, dimensionless (0 q 5 deg) 

1. 0-0.9 (0-26) 


0.9-0.8 (26-37) 


0.8-0.7 (37-46) 


0.7-0.6 (46-53) 


“Single-cap bin. 

b Assumed symmetry about principal plane. 


TABLE II. DATA SET RANGES AND POPULATION 


Desert 

site 

Spacecraft 

0 range, 
deg 

(j> range, 
deg 

0 O range, 
deg 

Number of 
measurements 

Sahara 

Nimbus 7 

6-87 

0-360 

0-53 

191 bins a 

Gibson 

Nimbus 7 

18-85 

5-336 

4-50 

63 points 

Saudi 

Nimbus 7 

6-86 

2-356 

4-47 

62 points 

Saudi 

Nimbus 6 

9-85 

1-358 

6-20 

51 points 


“Sufficient data were unavailable for 5 of 196 potential bins. 
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TABLE III. REGRESSION PARAMETERS FOR REFLECTANCE MODEL 


Desert site 

Spacecraft 

Vo 

Yi 

TV 

Csw 

a/Y 

Sahara 

Nimbus 7 

0.011 

0.920 

1.764 

0.33 

0.054 

Gibson 

Nimbus 7 

.009 

.623 

1.786 

.60 

.077 

Saudi 

Nimbus 7 

.008 

1.088 

1.678 

.18 

.057 

Saudi 

Nimbus 6 

.009 

1.186 

1.677 

.18 

.067 


TABLE IV. REGRESSION PARAMETERS FOR EMITTANCE MODEL 


Desert 

site 

Spacecraft 

Time 

Uo 

Clw 

M 

Llw, 0° . 

W/m 2 -sr 

Elw, 

W/m 2 

Sahara 

Nimbus 7 

11/78 

0.95 

0.01 

0.144 

113 

331 



to 

.85 

.01 

.117 

107 

317 



5/80 

.75 

.01 

.107 

101 

301 




.65 

.01 

.095 

95 

285 

Gibson 

Nimbus 7 

12/78 

0.99 

0.04 

0.170 

120 

348 



7/79 

.65 

.04 

.121 

98 

290 

Saudi 

Nimbus 7 

7/79 

0.99 

0.02 

0.164 

116 

337 



12/78 

.72 

.02 

.148 

104 

304 

Saudi 

Nimbus 6 

8/75 

0.98 

0.02 

0.176 

111 

321 
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